Electronic structure of the noncentrosymmetric 
superconductor MgioIrigBig 



B Wiendlocha, J Tobola and S Kaprzyk 

Faculty of Physics and Applied Computer Science, AGH University of Science 
and Technology, Al. Mickiewicza 30, 30-059 Cracow, Poland 

E-mail: bartekwOf atcat . ft j .agh.edu.pl 

Abstract. Electronic structure of a novel superconducting noncentrosymmetric 
compound MgioIrigBie was calculated using the Korringa-Kohn-Rostokcr 
method. Electronic part of the electron-phonon coupling constant, McMillan- 
Hopfield parameters, were calculated using the rigid-muffin-tin approximation 
(RMTA). The magnitude of the electron-phonon coupling constant A, analysing 
atomic contributions, is discussed. Our results show, that superconductivity in 
MgioIri9Bi6 is presumably mediated by electron-phonon interaction. 



PACS numbers: 74.25.Jb 

1. Introduction 

The interest in noncentrosymmetric structures, exhibiting superconductivity, has 
grown up in the past years. There are only a few examples of this type, which 
belong to various classes of materials, e.g. an antiferromagnetic heavy fermion 
system CePtaSi [I], ferromagnetic uranium compound Ulr, superconducting under 
pressure [2J, or non-magnetic ternary borides Li2Pd,3B [3] and L^PtaB [¥]. The main 
reason, why these systems are especially attracting, is related to the role of inversion 
symmetry in electron pairing. The absence of inversion symmetry may suppress the 
triplet pairing or mix singlet and triplet symmetry [3 13 13 • 

Very recently, Klimczuk and co-workers [5] synthesised a new type of intcrmetallic 
light-element based compound MgioIrigBi6, exhibiting superconductivity near 5 K. 
This novel material also belongs to the rare noncentrosymmetric structures, and 
crystallises in a large and rather complex bec cell (space group I-43m) [8]. MgioIrigBi6 
is partly similar to the Li2(Pd,Pt)3B system, since these structures contain an alkali 
metal (Li, Mg), boron, and heavy transition metal (Pd, Pt, Ir). 

In this work we intend to start the discussion on superconductivity mechanism in 
this unusual compound, analysing the electronic structure and the strength of electron- 
phonon coupling (EPC). Assuming the BCS-type behaviour in MgioIrigBi6, we study 
whether superconductivity is driven by light boron sublattices, like e.g. in MgB2, or 
by heavy transition metal atoms, as suggested for L^PdsB [5]. 

1.1. Computational details 

Electronic structure calculations were performed using the Korringa-Kohn-Rostoker 
(KKR) multiple scattering method [13 [TT] . The crystal potential was constructed in 
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the framework of the local density approximation (LDA) , using von Barth and Hedin 
formula |12] for the exchange-correlation part. For all atoms angular momentum 
cut-off Imax = 3 was set, k-point mesh in the irreducible part of Brillouin zone 
(IRBZ) contained over 200 points. Densities of states (DOS) were computed using the 
tetrahedron k-space integration technique, generating over 600 tetrahedrons in IRBZ. 
Due to the high atomic number of iridium (Z = 77) semi-relativistic calculations were 
performed, but neglecting the spin-orbit coupling, which is commented at the end of 
this paper. 




Figure 1. Unit cell of MgioIrigBig. The lack of inversion centre is clearly 
seen e.g. along the main diagonal, where Bl and Mgl atoms break the inversion 
symmetry. Generated by XCRYSDEN (27] , 

As far as the crystal structure is concerned, experimental lattice constant 
a = 10.568 A and atomic positions [5] were accounted for the computation (for clarity 
also shown in tabled]). Atoms in the unit cell were surrounded by muffin-tin (MT) 
spheres with following radii: Rm 9 = 2.82, Rj r = 2.50, Rb = 1-40 (in atomic units), 
filling about 60% of the cell volume. In the primitive cell of this system, 45 atoms 
occupy 7 inequivalent sites, which all are listed in tabled The noncentrosymmetricity 
of this system is important, and cannot be regarded as the effect of a lattice distortion, 
as observed e.g. in TJIr [5]. The inversion symmetry is broken by both boron sublattices 
(Bl and B2), as well as Mgl and Ir3 sites, thus the crystal has a half of symmetry 
operations of the cubic group Oh, i.e. only 24 operations. Among of all sublattices, 
the positions of iridium atoms are the closest to have inversion symmetry. Since Ir3 
occupies (x, z, z) sites, with x ~ 0.07, the full cubic symmetry is restored after shifting 
this position to (0, z, z). 

The analysis of superconducting properties is based on the computed McMillan- 
Hopfield (MH) parameters 77 [HI [14] which determine electronic part of the electron- 
phonon interaction, and directly enter the formula for the electron-phonon coupling 
constant A: 
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Figure 2. (a) Total DOS of MgioIrigBie . (b) Zoom near the Fermi level 
(Ep = 0) with atomic contributions. 



In equation ([T]), r\i is the MH parameter for each nonequivalent atom i which is 
characterised by the atomic mass Mj and averaged squared vibrational frequency (u)f). 
This equation divides the electron-phonon coupling constant into site-dependent parts, 
thus it allows to justify which sublattice gives the most important contribution to the 
total A. Calculations of MH parameters were preformed using the rigid muffin tin 
approximation (RMTA) [HI [TH1 E3 EH] , which gives the following expression for rj at 
each atomic site i: 

2 
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Here, I is the angular momentum number, n\(E F ) is the partial (angular-decomposed) 
density of states per spin at the Fermi energy (Ep), R\(r) is a regular solution of the 
radial Schrodinger equation, normalised to unity inside the MT sphere of radius R\ IT , 
and Vi(r) is the self-consistent, spherically-symmetric potential. N(Ep) is the total 
DOS at Ep per spin and per cell. For more detailed discussion of approximations 
involved in this method, see e.g. [THl HOj and references therein. 



2. Results and discussion 



The calculated total and atomic densities of states functions are presented in figure [2] 
As one could expect, electronic structure is dominated by iridium 5d states. The total 
DOS curve may be characterised as a collection of numerous van Hove singularities, 
reflecting the large number of atoms in the unit cell and various interatomic distances. 
The site-decomposed DOS are shown in figure [H and their values at the Fermi level 
are gathered in table[TJ The total DOS at Fermi level is about N(Ep) ~ 150 Ry -1 per 
formula unit. The densities at Ep, calculated per atom, are rather low (average value 
per atom: 3.3 Ry -1 ), and generally Ep is located outside of dominating peaks of DOS, 
on a small decreasing slope of Ir DOS. Strong hybridisation of d states of Ir and p 
states of B and Mg is manifested in the separation of bonding and anti-bonding states, 
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Table 1. Densities of states n(Ep), Z-decomposed DOS n; (Ry _1 /spin per atom) 
McMillan-Hopfield parameters w, MH parameters for each scattering channel 
Vl l+i (mRy/a.u. 2 per atom), and summary MH parameters for each site Wi X r\i 
(mRy/a.u. 2 per site). Wi is the number of atoms occupying site in primitive cell. 
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Wi 


n(E F ) 


n s 


iip 




Tlf 
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n sp 


v pd 


rf s 


Wi X f]i 


Irl 
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2.77 


0.15 


0.88 


1.86 
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2.8 


0.1 


1.4 


1.3 


2.8 


Ir2 
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3.16 


0.07 


0.34 


2.71 


0.03 


3.3 


0.0 


0.8 
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19.8 


Ir3 


12 


2.41 


0.06 


0.34 


1.97 


0.03 


2.6 


0.0 


0.6 


2.0 


31.2 


Bl 
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0.40 


0.03 


0.35 


0.01 


0.00 


1.2 


0.0 


1.2 


0.0 


4.8 


B2 


12 


0.27 


0.04 


0.22 


0.01 


0.00 


0.6 


0.0 


0.6 


0.0 


7.2 


Mgl 
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0.98 


0.17 


0.69 


0.10 


0.03 


0.7 


0.5 


0.2 


0.0 


2.8 


Mg2 
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0.09 


0.41 


0.14 


0.03 


0.2 


0.1 


0.1 


0.0 


1.2 



and Ep is located in the DOS valley (figure [2] and figure[3]). The interatomic distances, 
listed in table supports the enhanced p-d hybridisation, especially between Ir and 
B atoms (the smallest distances, 2.1 - 2.2 A). 

It is also interesting to compare the computed site-decomposed Ir densities, to 
the DOS of metallic fee iridium, which is presented in figure[3]Jd). The shape of lower 
part of DOS in Ir2 and Ir3 is quite similar to the case of Ir-fcc. This is probably due 
to the fact, that the Ir2-Ir3 and Ir3-Ir3 coordination, as well as interatomic distances, 
are very close to the fee phase (Ir-Ir distance is ~ 2.7 A in the aforementioned cases, 
see table [2]). However, the n(Ep) values on Ir atoms are much lower, comparing 
to fee structure (n(Ep) ~ 6.3 Ry _1 /spin in Ir-/cc), being the effect of enhanced 
hybridisation near Ep. Noteworthy, Irl has a quite different atomic coordination, 
with respect to Ir2 and Ir3 positions, being surrounded practically only by 4 boron 
atoms (Bl). It is clearly reflected by the apparently different DOS shape below Ep. 
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Figure 3. Site-decomposed DOS of Mgi Iri 9 Bi 6 (E F = 0). The total DOS of 
fee iridium metal is given for comparison in panel (d). 
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Table 2. The smallest interatomic distances between atoms (in A) and atomic 
positions in MgioIrigBig. The Ir-Ir distance in Ir-/ce is about 2.7 A, as between 
Ir2-Ir3 and Ir3-Ir3. 
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Irl 9.1 
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(2a) 
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(24g) 


Bl 






3.4 


3.3 


3.0 


3.0 


(0.3331, 0.3331, 0.3331) 
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2.6 
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Mgl 










5.0 


3.1 


(0.1127, 0.1127, 0.1127) 


(8c) 


Mg2 












3.1 


(0.1639, 0.1639, 0.4140) 


(12e) 



Electronic dispersion curves near the Fermi level are presented in figure [5] The 
bands located below —0.1 Ry (partly showed) are very flat, and form the narrow large 
DOS peaks, seen already in figure [U Conversely, the bands that cross Ep are quite 
dispersive, which results in the low DOS in this energy range. Noteworthy, there is an 
energy gap along the P-N direction (parallel to the k z axis in the reciprocal space). 
This may suggest some anisotropic transport properties of this compound. 

The calculated McMillan-Hopfield factors for all sites, with contributions from 
each scattering channel (I — ► l+l), are presented in table[TJ Estimation of the electron- 
phonon coupling constant A, using MH parameters, requires also the knowledge 
of average phonon frequencies. For such a large structure, containing 90 atoms 
in the cubic unit cell, phonon spectra calculations are difficult to be carried out. 
Nevertheless we can try to investigate the strength of electron-phonon interaction 
assuming reasonable values of (tuf) and studying their influence on estimated A values. 
Debye frequencies O of monoatomic crystals of iridium (fee), boron (rhomboedral) 
and magnesium (hep) may be helpful for choosing sensible phonon frequency range 
for our discussion. At first we assume, that the same type of atoms, occupying 
different sites in the cell, have similar average vibrational frequencies. Iridium, as 
the heaviest element, is expected to have the lowest (cof), Debye frequency of metallic 
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Figure 4. Dispersion curves near the Ep = in MgioIrigBig. 



Electronic structure of the noncentro symmetric superconductor MgioIrigBiQ 6 




15 20 25 30 35 40 40 50 60 70 80 90 100 15 20 25 30 35 40 

2 1/2 2 1/2 2 1/2 

<u h '> (meV) <w B > ~ (meV) <u Mg > (meV) 



Figure 5. Contributions to the electron-phonon coupling constant A in 
MgioIri9Big from iridium (a), boron (b) and magnesium (c) as a function of a root 
of average square frequency. Top curve on each graph is a sum of contributions 
from particular sites. Vertical lines mark the value of Si/V2. 



iridium is about 0/ r ~ 36 meV (420 K) [ST]. In contrast, the lightest boron is 
certainly expected to have the highest phonon frequencies, and for crystalline boron 
Qb — 100 meV (1200 K) [25] . Finally, magnesium average frequencies are expected 
to locate between the values of iridium and boron. Debye frequency of crystalline Mg 
is rather low: Qmq — 34 meV (400 K) J3JJ. We may also recall, that average square 
phonon frequency is often estimated as {uj 2 ) ~ which is a good approximation 

in monoatomic structures. In our case, these values may also be helpful for choosing 
feasible range of (u)f). 

Figure [5] presents estimated electron-phonon coupling constant A, associated with 
particular crystal sites, plotted as a function of average square phonon frequency. For 
each type of atoms, a wide frequency range was chosen, to illustrate the changeability 
of A. In the case of iridium, the largest contribution to A comes from Ir3 sublattice, 
due to the large population of this site. Among boron and magnesium sites, the B2 
and Mgl atoms provide the largest contributions. However, one has to remember that 
this comparison is valid only if we assume identical frequencies for the same atoms 
at different sites. In this simplified analysis, we may also plot the overall coupling 
constant for constituent atoms, by adding the contributions from each site, which is 
illustrated by solid lines in figure [5] 

As we can see, the obtained partial coupling constants are not high for each atom. 
At the moment, we are not aware of any experimental findings of the EPC constant 
in this compound. However, we can try to estimate the range of "experimental" A, 
analysing the magnitude of the observed critical temperature. If we assume, that 
we are dealing with BCS-type superconductivity, we may substitute the experimental 
value of T c = 4.5 K into the McMillan formula for T c [15] : 

1.04(1 + A) 



rp e 

Tc = I¥ exp 



(3) 



A- M *(1 + 0.62A) 

However, because the Debye temperature of MgioIrigBi6 is not known yet, we 
plot A in figure E] as a function of 0, for typical values of Coulomb pscudopotential 
parameter /j*. The resulting EPC constant A varies between 0.5 and 0.75, and for 
e.g. = 250 K we obtain A ~ 0.60, whereas for = 350 K we get A ~ 0.55 
(for jj* = 0.11), as one can see in the figure [S] Thus, if MgioIrigBi6 is treated as 
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Figure 6. "Experimental" value of the electron-phonon coupling constant A, 
evaluated from the McMillan formula, using observed Tq = 4.5 K, plotted as a 
function of Debye temperature © for three values of Coulomb pseudopotential 
parameters fi* . 



a conventional superconductor, the lower limit of A is expected to be 0.5. In order 
to get this value from our RMTA calculations, relatively low average frequencies for 
all constituent atoms are required. If we take the following \J (cuf } : 14 meV for Ir, 
50 meV for B and 20 meV for Mg, we obtain A/ r ~ 0.3, As ~ 0.1 and \m§ — 0.1, 
which gives the expected lower limit of EPC constant (A ~ 0.5). If higher values of A 
are experimentally observed, it will indicate cither additional contributions, neglected 
within the RMTA framework, or even lower phonon frequencies, than the values used 
in the aforementioned estimations. Especially iridium contribution is sensitive to the 
change of \J {oSj r ) parameter, i.e. increasing it to 25 meV (probably an upper limit) 
results in decrease of A/ r to 0.1, that is to the value found for other sublattices. 
Generally, our calculations indicates, that iridium sublattices seem to be the most 
important for the onset of superconductivity in MgioIrigBig. 

Finally, we shortly comment the possible influence of the lack of inversion 
symmetry and spin-orbit (SO) interaction on the electronic structure and 
superconductivity. The absence of a centre of inversion in a superconductor affects 
the symmetry of superconducting state, allowing for an admixture of singlet and 
triplet components [HI [7] . Because the singlet pairing is based on the time- reversal 
symmetry |23j . which is present as long as the compound is not magnetic, mainly 
the triplet channel is affected by the lack of inversion, and the superconductivity may 
be even suppressed, see e.g. [5j El [7]. The SO interaction, which mixes the initial 
spin-up and spin-down electronic states, plays an additional role. It was found, that 
it may control the mixing of parity of the superconducting state [7] . This seems to be 
the case of I^PdaB and L^PtaB compounds, where specific heat [H] and NMR [35] 
measurements strongly support phonon-mediated isotropic superconductivity, while 
penetration depth measurements suggests an admixture of spin-singlet and triplet 
components in the superconducting energy gap [26] , with larger triplet component in 
the Pt-case. This kind of experimental study for MgioIrigBi6 should be prior to the 
theoretical discussion of gap symmetry in this compound. 

As far as the band structure of MgioIrigBi6 is concerned, the modifications due 
to the SO interaction are not expected to significantly affect the obtained values of 
DOS and MH parameters. For metallic fee iridium, our value of MH parameter, 
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calculated also neglecting SO interaction, rjj r ~ 135 mRy/(a.u.) 2 , gives the correct 
magnitude of the electron-phonon coupling constant A = 0.32, comparing to the 
observed A = 0.34 13J . Here, the average square frequency was estimated from the 
formula (w 2 ) ~ ^6 2 , using the Debye temperature = 420 K. 

3. Summary and conclusions 

The results of LDA electronic structure calculations of new MgioIrigBi6 
superconductor were presented. The main contributions to densities of states 
near Ep are provided by iridium atoms. The electron-phonon coupling constant 
A was roughly estimated, using the calculated McMillan-Hopficld parameters and 
qualitative discussion of average phonon frequencies. We discussed the relation of 
the experimental transition temperature and the magnitude of EPC coupling. Within 
the rigid-muffin-tin approximation, the main contribution to A comes from iridium, 
with smaller contributions from boron and magnesium. If any information about 
the dynamic properties of atoms in this compound become available, our analysis 
presented of figure will allow to find better theoretical estimation of A. The location 
of Ep on the slope of Ir DOS peak leads to rough suggestion, that hole doping on 
iridium sites, e.g. with rhodium, may increase the densities and MH parameters. 
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